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Abstract 

Purpose When performing uncertainty propagation, most 
LCA practitioners choose to represent uncertainties by sin¬ 
gle probability distributions and to propagate them using 
stochastic methods. However, the selection of single proba¬ 
bility distributions appears often arbitrary when faced with 
scarce information or expert judgement (epistemic uncer¬ 
tainty). The possibility theory has been developed over the 
last decades to address this problem. The objective of this 
study is to present a methodology that combines probability 
and possibility theories to represent stochastic and epistemic 
uncertainties in a consistent manner and apply it to LCA. A 
case study is used to show the uncertainty propagation 
performed with the proposed method and compare it to 
propagation performed using probability and possibility 
theories alone. 

Methods Basic knowledge on the probability theory is first 
recalled, followed by a detailed description of epistemic un¬ 
certainty representation using fuzzy intervals. The propaga¬ 
tion methods used are the Monte Carlo analysis for probability 
distribution and an optimisation on alpha-cuts for fuzzy in¬ 
tervals. The proposed method (noted as Independent Random 
Set, IRS) generalizes the process of random sampling to 
probability distributions as well as fuzzy intervals, thus mak¬ 
ing the simultaneous use of both representations possible. 
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Results and discussion The results highlight the fundamental 
difference between the probabilistic and possibilistic represen¬ 
tations: while the Monte Carlo analysis generates a single 
probability distribution, the IRS method yields a family of 
probability distributions bounded by an upper and a lower 
distribution. The distance between these two bounds is the 
consequence of the incomplete character of information 
pertaining to certain parameters. In a real situation, an excessive 
distance between these two bounds might motivate the 
decision-maker to increase the infonnation base regarding cer¬ 
tain critical parameters, in order to reduce the uncertainty. Such 
a decision could not ensue from a purely probabilistic calcula¬ 
tion based on subjective (postulated) distributions (despite lack 
of information), because there is no way of distinguishing, in 
the variability of the calculated result, what comes from true 
randomness and what comes from incomplete information. 
Conclusions The method presented offers the advantage of 
putting the focus on the information rather than deciding a 
priori of how to represent it. If the information is rich, then a 
purely statistical representation mode is adequate, but if the 
information is scarce, then it may be better conveyed by 
possibility distributions. 

Keywords Confidence index ■ Distribution ■ Fuzzy sets • 
Intervals ■ Possibility ■ Probability • Uncertainty 
propagation ■ Uncertainty representation 

1 Introduction 

Life cycle assessment (LCA) aims at modelling complex 
systems that usually encompass a number of compartments 
of the biosphere and the technosphere. Results rely on 
several choices and large amounts of data are affected by 
uncertainty. These uncertainties have been described exten¬ 
sively, e.g., by Reap et al. (2008) and Williams et al. (2009). 
Characterising and assessing uncertainties is important to 
make decision support models more transparent, robust and 
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reliable. Uncertainty analysis gathers numerous methods 
with different means and goals from qualitative assessment 
to sensitivity analysis and uncertainty propagation; see 
Morgan and Henrion (1990) for an overview on uncertainty 
analysis and Clavreul et al. (2012) for a tiered approach to 
uncertainty analysis in LCA applied to waste management. 

The focus of the present study is on uncertainty propa¬ 
gation which aims at quantifying the uncertainty of the 
results of an LCA study. Uncertainty propagation can be 
performed using different uncertainty representations and 
propagation methods. With respect to parameter uncertainty, 
the common practice in LCA consists in representing un¬ 
certain parameters by single probability distributions, e.g. a 
nonnal distribution characterized by an average and a stan¬ 
dard deviation. Databases such as the ecoinvent database 
(Frischknecht et al. 2005) rely increasingly on probability 
distributions to represent parameter uncertainty. The most 
commonly used method to propagate probability distribu¬ 
tions is the Monte Carlo analysis, as shown by Lloyd and 
Ries (2007) who reviewed 24 LCA studies that included 
uncertainty analysis. This method is implemented in many 
calculation tools and consists in randomly sampling values 
in the probability distributions of input parameters, to obtain 
the frequency distribution of the calculated results. 

Flowever, a fundamental problem of the probabilistic repre¬ 
sentation lies in the selection of probability distributions when 
faced with scarce information or expert judgement. The review 
by Lloyd and Ries (2007) showed that the choice of probability 
distributions is often poorly justified and relying on estima¬ 
tions. Yet, the result of the uncertainty propagation is totally 
depending on the ‘a priori’ defined probability distributions. 
Bayesian methods (Lindley 1971) could address this shortcom¬ 
ing by updating these prior distributions based on new data and 
Bayes’ theorem of conditional probabilities. Flowever, this is 
almost never implemented in LCA due to the impossibility of 
measuring and validating results. This introduces confusion 
between the two distinct natures of uncertainty: truly stochastic 
uncertainty which refers to variability of data e.g. in time, space 
and technology, and epistemic uncertainty related to our lack of 
knowledge, e.g., due to measurement errors or to an insufficient 
number of measurements. While classical probability theory 
was developed to address stochastic uncertainty (i.e. related to 
variability and fluctuations), more recent information theories 
are required to address incomplete/imprecise information 
(Dubois and Prade 2009). 

To address this problem and handle modelling in presence 
of imprecise infonnation, possibility theory has been developed 
over the last decades (Dubois and Prade 1988). The simplest 
example is the representation of parameters as min-max in¬ 
tervals instead of crisp (precise) numbers, as used by Chevalier 
and Le Teno (1996). The concept can be extended to fuzzy 
intervals (possibility distributions) which express preferences 
within intervals. More detailed presentation of fuzzy intervals is 


provided in the “Methods” section. Fuzzy intervals have been 
first applied to the field of LCA to save time and costs by 
avoiding the need for precise quantification of flows, e.g., by 
Weckenmann and Schwan (2001) and Gonzalez et al. (2002). 
Fuzzy linguistic descriptors have also been used to calculate 
life cycle inventories (LCI) and evaluate data quality (Ardente 
et al. 2004), to normalise and weigh characterised impacts 
(Guereca et al. 2007) and to support interpretation of results 
and ranking alternatives using multi-criteria analysis (Benetto 
et al. 2008). Besides, in an LCA model dedicated to fuel 
evaluation (namely POLCAGE), Tan et al. (2004) represented 
parameter uncertainties using possibility distributions and 
propagated them using fuzzy arithmetics. Tan (2008) 
formalised the integration of fuzzy intervals into a matrix- 
based LCI model, supported later by a mathematical proof by 
Heijungs and Tan (2010) re-examined by Cruze et al. (2013). 
Finally, Andre and Lopes (2012) proposed to enhance the 
mathematical and physical understanding of the application 
of possibility theory to LCA, by providing clear definitions of 
terms and comparing the possibility and probability represen¬ 
tations and propagation results. 

The objective of this study is to present a methodology 
that combines probability and possibility theories to repre¬ 
sent stochastic and epistemic uncertainty in a consistent 
manner, and apply it to LCA. The method is compared to 
uncertainty propagation performed with probability and 
possibility theories alone, using a case study where global 
warming benefits associated with bioenergy from energy 
crops cultivation are assessed. 

2 Methods 

This section describes how probability and possibility theo¬ 
ries can be used to represent uncertainties and propagate 
them through a model. A joint-propagation method, pro¬ 
posed by Baudrit et al. (2006), is presented and applied to 
a case study. In addition, two other propagation methods are 
also applied for the purpose of comparison. 

2.1 Probability theory 

A probability is a measure of the likelihood that an event will 
occur. A probability distribution describes the probabilities of 
different outcomes of a statistical experience: for a random 
variable X, a probability distribution gives for each value x the 
probability P(x) that X takes the value x. Probabilities follow 
certain rules: they take only values between 0 and 1 and the 
sum of the probabilities of all possible outcomes is 1. In the 
case of continuous variables, probability distributions are 
often represented by their cumulative distribution function 
(cdf), the probability that X be less than x: F{X)~P\X<x\ 
Another representation, the probability density function 
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(noted here as pdf), can be obtained by deriving the cdf. It 
represents the density of probabilities: for some small incre¬ 
ment Ax, fix). Ax is the probability that A'falls in the interval of 
length Ax around x (Morgan and Henrion 1990). In theory, 
selection of a probability distribution should be based on a 
sufficient amount of data to allow a statistically representative 
assessment of the parameter’s variability. However, in the 
context of LCA, this is often not technically feasible and the 
selection of a distribution often relies on partial information 
(scarce measurements) or on expert judgement. 

Uncertainty propagation of probability distributions can be 
perfonned by different methods, the most common one being 
the Monte Carlo analysis, as used by, e.g. Huijbregts et al. 
(2003) or Sonnemann et al. (2003). In this analysis, a value is 
randomly sampled for each parameter in its distribution and 
by using the obtained set of values, the model result is calcu¬ 
lated. By repeating this operation a sufficient number of times, 
a cdf is obtained for the result. Other sampling methods are 
more adapted to large data sets, such as Latin hypercube, as 
used by Thabrew et al. (2008). Finally, calculations can also 
be performed analytically using Taylor series expansions to 
approximate the result’s uncertainty, as implemented by Hong 
et al. (2010) and Imbeault-Tetreault et al. (2013). 

2.2 Epistemic uncertainty representation 

As shown by several authors (e.g. Ferson and Ginzburg 1996), 
there is a fundamental difference between tale random vari¬ 
ability, as depicted by a single probability distribution, and 
epistemic uncertainty, due to incomplete or imprecise informa¬ 
tion. Possibility theory (e.g. Dubois and Prade 2008) provides 
a framework to address this type of information. Possibility 
theory assigns degrees of likelihood (possibility) to intervals of 
values rather than precise values, yielding a fuzzy interval (or 
fuzzy number or fuzzy set). The simplest fuzzy interval is the 
well-known min-max interval. If the parameters involved in a 
model are represented by intervals, interval propagation can be 
performed using either interval calculus in the case of simple 
models or else an optimisation algorithm in the case of more 
complex models. In a Bayesian framework, application of the 
principle of maximum entropy to interval-type information 
leads to selecting a uniform probability distribution between 
the limits of a min-max interval (e.g. Shulman and Feder 
2004). But this results in selecting only one distribution 
amongst all the possible probability distributions bounded by 
the following two cumulative distributions: 

PfiX) = 0 if X < min and P U (X) = 1 otherwise, 

P\(X) = 0 if X < max and Pi(X) = 1 otherwise. 

Where P u and P\ are the upper and lower limits of the 
family of probability distributions defined by the min-max 


interval. Selecting just one representative of the family of 
probability distributions introduces a bias in the analysis and 
a confusion between true variability (as depicted by a single 
distribution) and imprecision (as depicted by an interval). 

When richer information is available, the concept of in¬ 
tervals can be extended to fuzzy intervals (also called pos¬ 
sibility distributions) where preference is given to certain 
values (see Dubois and Prade 1988). In a possibility distri¬ 
bution, degrees of likelihood (possibility) between 0 and 1 
are assigned to specific parameter intervals. In the example 
depicted in Fig. 1, the most likely interval (the “core” of the 
possibility distribution), i.e. values between 18 and 20, is 
assigned a likelihood of unity, while values located outside 
the “support” of the distribution (i.e. values between 14 and 
23) are assigned a possibility of zero. Intervals selected at 
different levels of possibility, called alpha-cuts, correspond 
to confidence intervals with confidence 1-alpha. Thus, a 
possibility distribution yields a lower bound (P a > 1-ct) on 
the probability that the parameter value should lie within a 
given alpha-cut. As in the case of the simple min-max 
interval, a fuzzy interval can be depicted as a family of 
probability distributions, limited by an upper and a lower 
cdf, as shown in Fig. 1. While the function presented in 
Fig. 1 is a trapezoidal distribution, more complex distribu¬ 
tions can be adopted to suit available information (see 
Dubois and Prade 1988). Fuzzy intervals are particularly 
well suited for representing subjective judgements, com¬ 
monly used in most LCA studies, because they adopt the 
language of experts, when describing the possible values a 
parameter can take in presence of incomplete information 
(Dubois 2006). If an expert is able to answer the following 
two questions: (1) can you provide a range within which you 
are confident that the parameter value should lie? and (2) 
can you express a preferred value or interval of values 



.Upper probability bound 

---- Lower probability bound 

Fig. 1 Example of a possibility distribution 
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within this range? Then the provided information can be 
formalized as a possibility distribution. 

2.3 Propagation methods 

The method used to propagate fuzzy intervals in the general 
case is very analogous to the Monte Carlo method using single 
probability distributions, except that in the case of parameters 
represented by fuzzy intervals, intervals are randomly sam¬ 
pled instead of single values, based on a-cuts. As shown 
above, for a given possibility distribution, an alpha-cut is an 
interval containing all values with a degree of possibility 
higher than alpha (0<a< 1). An example of o-cut is presented 
in Fig. 1 : for alpha=0.6 the a-cut is the interval [16.4; 20.9], If 
the model is relatively simple and monotonous, propagation 
of the fuzzy intervals through the model can be performed 
simply using interval calculus on alpha-cuts. For alpha=0 to 1 
with, e.g. step=0.1, the min and max values of the model are 
determined for all values of the alpha-cuts. However, if the 
model is not monotonous, it may not be possible to determine 
the min and max values of the model based solely on the min 
and max values of the alpha-cuts. In this case it is necessary to 
use an optimization algorithm to find the min and max values 
of the model for all parameter values within the alpha-cuts. 

If certain parameters are represented by fuzzy intervals 
while others are represented as single cdfs, the Monte Carlo 
method can be used to randomly sample the cdfs, while 
optimization on the alpha-cuts is performed in a second step 
(see Guyonnet et al. 2003 and Baudrit et al. 2005). Baudrit 
et al. (2006) developed a slightly different method (dubbed 
as the IRS, Independent Random Set method), whereby 
random sampling is performed on both the cdfs and the 
fuzzy intervals. Couso et al. (2000) showed that the IRS 
method is a systematically conservative counterpart to the 
calculation with random quantities under stochastic inde¬ 
pendence (classical Monte Carlo method on cdfs). The 
schematic of the IRS method, used herein, is the following. 

Given an LCA model with n parameters represented by 
probability distributions and m parameters represented by 
fuzzy intervals, 

1. Generate n+m random numbers between 0 and l:x\, x 2 , 

• • 'X n + m . 

2. Sample the n probability distributions to obtain n ran¬ 
dom variables: p\,p 2 ... p„- 

3. Sample the m fuzzy intervals to obtain m intervals: /), 

I 2 ’ • • 

4. Calculate the smallest (Inf) and largest (Sup) values of 
the LCA result obtained for all combinations of values 
contained in the m intervals I,. 

5. Return to step 1 and repeat to times. 

6. Obtain the probability bounds of the LCA results from 
the co Inf and Sup values as shown below. 


The IRS method yields a random interval made up of u> 
intervals. This random interval is then summarized in the form 
of a pair of upper and lower cdf (see Baudrit et al. 2005) using 
the Plausibility and Belief functions of the theory of evidence 
(Shafer 1976). This theory assigns probability weights (noted 
m) to intervals (called focal sets; A;) instead of simply point 
values (the limiting case of a classical probability distribu¬ 
tion). Considering the proposal (noted B) “LCA result lies 
below a specified target level”, the probability that this pro¬ 
posal is true is comprised between the degree of Plausibility 
(an upper bound on probability) and the degree of Belief (a 
lower bound on probability) defined by Shafer (1976) as: 

Bel ( s > = E, (i) 

«( s ) = E, w .”W) P) 

Bel(R) is thus the sum of the weights of all subsets A, (/= 
1 to n where n is the number of subsets) such that A t is 
completely included within prescribed set B, while Pl(Zi) is 
the sum of the weights of all subsets such that the 
intersection of A, and B is non empty. In other words, 
Belffi) gathers the imprecise evidence that asserts B, while 
P1(R) gathers the imprecise evidence that does not contradict 
B. The interval [Belffi), P1(R)] contains all potential proba¬ 
bility values induced by the mass function m. In practice, PI 
is obtained by ordering the u> Inf values in increasing order, 
and assigning a frequency 1 ho to each value, while Bel is 
obtained likewise from the Sup values. These functions will 
be depicted graphically in the application section below. 

2.4 Interpretation of results in a decision-making framework 

If at least one parameter in a model is represented by a fuzzy 
interval, the uncertainty propagation will result in a family 
of probability distributions (delimited by the PI and Bel 
functions described previously), rather than in a unique 
probability distribution. As suggested by Dubois and 
Guyonnet (2011), this may prove impractical in a decision¬ 
making framework. These authors therefore propose to 
compute a single distribution as a weighted average of the 
upper and lower distributions, with the selected weight 
reflecting the decision-makers attitude with respect to risk. 
The resulting distribution, referred to as a “confidence in¬ 
dex” by Dubois and Guyonnet (2011), is computed from: 

/(«;> b t ) = la, + (7 -l) A (3) 

where a, and b, are the limits of the interval defined at 
probability level i. This approach, which is based on earlier 
work by Hurwicz (1951), thus computes a single indicator 
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as a weighted average of focal element bounds. It 
achieves a trade-off between upper and lower probability 
bounds which, following the context, will constitute ei¬ 
ther optimistic or pessimistic estimates. While it is rec¬ 
ognized that the choice of weight a is subjective, it 
should be underlined that this subjectivity is only intro¬ 
duced at the decision-making step in the form of a single 
cdf used as a sensible reference displayed along with the 
pessimistic and optimistic outputs. This approach is very 
different from displaying a single distribution obtained by 
propagating single distributions selected arbitrarily at the 
beginning of the risk analysis step. 

3 Case study 

3.1 Goal and scope 

The objective of this LCA study is to exemplify and apply 
uncertainty propagation using different hypotheses with re¬ 
spect to input parameter uncertainty, for the purpose of com¬ 
parison. To this end, a specific LCA case study was selected in 
order to show the differences between the probability and 
possibility theories and how they can be combined in order 
to better represent uncertainties in LCA. The selected case 
study investigated the environmental sustainability of willow 
cultivation for bioenergy production through co-firing in 
large-scale combined heat and power (CHP) plants, based 
on results from Tonini et al. (2012) (see “Background—case 
study”). Emphasis was placed on how the different uncer¬ 
tainties associated with the inventory data can be represented 
based on available knowledge (e.g., from measurements, lit¬ 
erature or expert estimates). The uncertainty associated with 
the environmental impact of the system was quantified with 
each individual uncertainty method to identify the major dif¬ 
ferences between them and recommend a best practice. The 
focus of this study was on the global warming impact catego¬ 
ry. The functional unit was the cultivation of 1 ha of Danish 
land for bioenergy production (CHP). The geographical scope 


was Denmark and the temporal scope 20 years. Figure 2 pre¬ 
sents the processes included in the LCA system boundary. 

3.2 Background—case study 

Reduction of fossil fuel consumption in the energy sector 
through increase of fluctuating renewable energy sources (e.g. 
wind energy and hydropower) and bioenergy is a fundamental 
step towards more sustainable energy systems (Tonini and 
Astrup 2012). However, biomass resources available for 
bioenergy are limited as biomass is already used today for a 
number of purposes (e.g. animal feeding and bedding, improve¬ 
ments of agricultural soil, etc.). Thus, cultivation of energy 
crops for bioenergy production may be needed. One of the 
most critical impacts associated with energy crops is related 
to land use changes (LUC) defined as the consequences deter¬ 
mined by the conversion of the land from one use to another use 
(Edwards et al. 2010; Searchinger et al. 2008; Searchinger 
2010). LUC are distinguished between direct (dLUC) and 
indirect (iLUC). The dLUC impacts are associated with the 
consequences of cultivating the selected energy crops in place 
of an established food crop. The iLUC impacts are related to the 
consequences of converting land presently not used for crop 
cultivation to cropland, as a result of the induced demand for 
the initially displaced food crop. In order to evaluate the envi¬ 
ronmental sustainability of bioenergy systems, LCA is often 
used. For instance, in Tonini et al. (2012), a case study based on 
cultivation of three perennial crops (ryegrass, willow and 
Miscanthus ) in Denmark was presented. The authors compared 
the environmental performance of anaerobic digestion, gasifi¬ 
cation, direct combustion and co-firing. With respect to global 
warming, co-firing of willow appeared to be the most environ¬ 
mentally sound option, though CCL savings were generally low 
as a result of LUC. 

3.3 Modelling and data 

The modelling of the bioenergy system (primarily CCB and 
N 2 0 flows) was based on the inventory data provided by 


Fig. 2 System boundary of the 
selected LCA case study 
(dashed lines : avoided 
processes) 
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recent studies: Hamelin et al. 2012 (cultivation of willow and 
of the marginal crop displaced, i.e. spring barley) and Tonini et 
al. 2012 (storage, pre-treatments, energy conversion processes 
and estimates of iLUC). Spring barley was assumed as the 
marginal crop, i.e. the food crop which would likely react to 
changes in demand or supply of energy crops (Dalgaard et al. 
2008; Schmidt 2008; Weidema 2003). Coal-fired power plants 
and natural gas-fired power plants were assumed as the mar¬ 
ginal technologies for respectively electricity and heat produc¬ 
tion (Energistyrelsen 2011; Weidema et al. 1999; Weidema 
2003). The overall environmental impact on global warming 
was thus calculated as the sum of the following processes (see 
Appendix for further details): 

I. Cultivation of willow; 

II. dLUC, i.e. the impacts/savings associated with the 
replacement of the marginal crop; 

III. iLUC, estimated after Tonini et al. (2012); 

IV. Storage, pre-treatments: emission of carbon dioxide; 


V. Co-firing: emissions of carbon dioxide; 

VI. Avoided energy production (i.e. avoided emissions of 
GHGs from fossil fuel combustion). 

It was assumed that all carbon degradation during drying, 
storage and combustion was in the form of carbon dioxide 
(methane emissions were negligible) and machinery-related 
processes such as fertilizer spreading and tillage were not 
included because they contribute to the results only to a 
minor extent (Tonini et al. 2012). Selected modelling data 
(referred to below as ‘parameters’) related to CO 2 and N 2 O 
flows throughout the bioenergy system were associated with 
uncertainty (after Hamelin et al. 2012 and Tonini et al. 
2012). The uncertainty representation modes are shown in 
Table 1. For the purpose of comparison, trapezoidal distri¬ 
butions were selected for both probability density function 
and fuzzy interval modes of representations. Their supports 
and cores are respectively delimited by the values [a, d\ and 
[b, c ] presented in Table 1 , estimated based on the different 


Table 1 Assumed parameter distributions, values rounded to two significant digits 


Description 

Unit 

Limits of the trapezoidal 

Source of information 

Preferred 



distributions 


representation 




a 

b 

c 

d 



Net carbon uptake from atmosphere, 

Mg C ha 1 year 1 

3 

6 

6 

9 

Tonini et al. 2012 

Fuzzy 

willow cultivation 

Net carbon uptake from atmosphere, 

Mg C ha 1 year 1 

1 

2 

2 

3 

Tonini et al. 2012 

Fuzzy 

barley cultivation 

N 2 O direct emissions, willow 

kg N ha 1 year 1 

0.8 

1.69 

1.69 

2.5 

Tonini et al. 2012 

Fuzzy 

cultivation 

N 2 0 direct emissions, barley 

kg N ha 1 year 1 

0.9 

1.9 

1.9 

2.9 

Tonini et al. 2012 

Fuzzy 

cultivation 

N 2 0 indirect emissions, willow 

kg N ha 1 year 1 

0.1 

0.22 

0.22 

0.33 

Tonini et al. 2012 

Fuzzy 

cultivation 

N 2 0 indirect emissions, barley 

kg N ha 1 year 1 

0.3 

0.56 

0.56 

0.8 

Tonini et al. 2012 

Fuzzy 

cultivation 

Indirect land use change 

Mg C0 2 -eq ha 1 

189 

398 

398 

610 

After Tonini et al. 2012 a 

Fuzzy 

Yield of cultivation of willow 

Mg DM ha~' 

8.7 

12.7 

12.7 

16.7 

Tonini et al. 2012 

Probability 

Yield of cultivation of barley 

Mg DM ha~‘ 

3.35 

4.85 

4.85 

6.35 

Tonini et al. 2012 

Probability 

Carbon content of willow 

% DM 

0.47 

0.48 

0.49 

0.50 

Tonini et al. 2012 b 

Probability 

Loss of carbon during storage 

% 

0.035 

0.048 

0.048 

0.061 

Tonini et al. 2012 

Fuzzy 

Lower heating value of dry 

GJ Mg~‘ DM 

16.7 

17.6 

19 

19.8 

Tonini et al. 2012 b 

Probability 

matter (willow) 

Water content of willow after 

% 

0.15 

0.2 

0.3 

0.35 

Tonini et al. 2012 

Fuzzy 

field drying 

Electricity recovery from LHV 

% 

0.35 

0.38 

0.38 

0.41 

Danish Energy Agency 
and energinet.dk (2010) 

Probability 

Heat recovery from LHV 

% 

0.44 

0.52 

0.52 

0.6 

Danish Energy Agency 
and energinet.dk (2010) 

Probability 

GHG emissions from electricity 
production in DK 

Mg C0 2 -eq MW1T 1 

0.66 

0.92 

0.92 

1.05 

Personal communication, DONG 
Energy A/S et al. (2010) 

Probability 

GHG emissions from heat 
production in DK 

Mg C0 2 -eq GJ -1 

0.04 

0.05 

0.06 

0.07 

Estimations based on the 
ecoinvent database 

Fuzzy 


a This includes the conversion of land and the effects of cultivating the reacting crop on newly converted land 
b Nineteen values extracted from articles and the Phyllis and Biodat databases referenced in Tonini et al. 2012 
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sources presented. The choice between probability and 
fuzzy interval representations was made based on the 
quantity and quality of data available for each individ¬ 
ual parameter. For example, Fig. 3 shows the 19 values 
collected in literature and databases for the lower 
heating value (LFIV) of willow. This significant amount 
of data enabled to define a trapezoidal distribution and 
to select a representation with probability distributions 
in the joint-propagation method. Conversely, very scarce 
information could be found on iLUC; therefore its un¬ 
certainty distribution was defined based on expert 
judgement and the representation using a possibility 
distribution was preferred in the joint-propagation 
method. 

Further, it is necessary to fix correlations when they are 
known, in order to avoid non-physical combinations of 
parameter values during random sampling process. 
Correlations between the following parameters were identi¬ 
fied and implemented: 

- Cultivation yield and net carbon uptake for willow, 

- Cultivation yield and net carbon uptake for barley, 

- Cultivation yields for willow and barley (as they depend 
highly on soil and climate properties), 

- LHV and carbon content of willow. 

These correlations were implemented by direct linear 
correlation: the carbon content was implemented as a func¬ 
tion of LHV while net carbon uptakes for willow and barley 
and the yield of barley cultivation were all implemented as 
functions of the yield of willow cultivation. Note that a 
fuzzy correlation could also be implemented, whereby the 
selection of one parameter generates an interval for the 
correlated parameter rather than a precise value (see 
Guyonnet et al. 2003). 

Heat and electricity recovery are assumed to be inde¬ 
pendent because the power plants are considered as 
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Fig. 3 Data collected in 19 studies for LHV of willow 


extraction condensing power plants. N 2 0 emissions (both 
direct and indirect) were considered as independent for 
willow and spring barley because they are linked to 
fertilizer use. 


4 Results 

Uncertainty propagation was performed using the fol¬ 
lowing three methods: Monte Carlo with cdfs, fuzzy 
calculation and the IRS method. The algorithm and its 
implementation in MATLAB 7 are provided as 
Electronic Supplementary Material. In the second and 
third uncertainty method, minimum and maximum 
values were calculated using the global search function 
of MATLAB 7. The cumulative distribution functions of 
the calculated results are presented in Fig. 4. The x-axis 
shows the impact of the system on global warming: a 
positive result means that the cultivation and co¬ 
combustion of willow contributed to more greenhouse 
gases emissions than current practice. 

The distribution obtained with Monte Carlo simulation 
(in full line) suggests a 65 % probability that willow culti¬ 
vation and combustion was beneficial compared to current 
practice. According to this simulation, the average benefit 
was -48 Mg C0 2 -eq ha” 1 for 20 years, with a standard 
deviation of 124 Mg C0 2 -eq ha -1 , and a 95 % confidence 
interval between -305 and 194 Mg C0 2 -eq ha” 1 . 

When implementing the same distributions as fuzzy in¬ 
tervals instead of probability distributions, two curves were 
obtained: a plausibility and a belief distribution (see Fig. 4). 
They are the respective upper and lower limits of the family 
of probability distributions obtained with fuzzy intervals. In 
this case study, the proposal evaluated was “The impact of 
willow cultivation for bioenergy production on global 
warming is below a specific target”. Thus the plausibility 
distribution represents here the most “optimistic” probabil¬ 
ity distribution: it is obtained from the most favourable 
values of input possibility distributions. On the other hand, 
the belief distribution represents the most “pessimistic” out¬ 
come achievable: the impact on global warming cannot be 
larger than this distribution, considering the input informa¬ 
tion. The global warming potentials resulting from the fuzzy 
calculus were between -603 and 400 Mg C0 2 -eq ha” 1 
(95 % confidence interval) and most likely between -240 
and 100 Mg C0 2 -eq ha” 1 . These very wide ranges result 
from the fact that the rich information available for some 
parameters was only modelled as fuzzy information in this 
calculation. 

In the third method, either mode of uncertainty represen¬ 
tation was selected, based on available information. Two 
distributions were again obtained, thus defining a family of 
distributions which again encompasses the purely 
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Fig. 4 Cumulative distribution 
functions of greenhouse gas 
emissions of cultivation and co¬ 
combustion of willow (in Mg 
CC> 2 -eq ha~’) obtained with 
three uncertainty propagation 
methods: Monte Carlo, fuzzy 
calculus and IRS method (1,000 
runs) 
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probabilistic result. Note that the distance between the upper 
and lower probability bounds, which directly reflects the 
incomplete nature of information regarding certain parame¬ 
ters, is less than in the case of the purely possibilistic 
calculation, because in this case certain parameters are rep¬ 
resented by single cdfs. 

Also depicted in Fig. 4 is the confidence index calculated 
by assigning a weight of 1/3 to the “optimistic” IRS result 
and 2/3 to the “pessimistic” result. Putting all the weight on 
the pessimistic bound would seem exaggerated, as it would 
be neglecting all information suggesting a more favourable 
outcome, while putting all the weight on the optimistic 
bound would appear as unrealistically biased. The selected 
weights of 1/3 and 2/3 are proposed as a “reasonably con¬ 
servative” compromise. 

5 Discussion 

In this study, it was recognized that the level of infor¬ 
mation was low for 10 out of the 17 parameters (cf. 
Table 1). The results highlighted the fundamental differ¬ 
ence between the probabilistic and possibilistic represen¬ 
tations: while the Monte Carlo analysis produces a crisp 
(precise) result on the probability of exceeding the 
baseline emissions (represented by a global warming 
potential of zero in Fig. 4), the IRS method yields a 
family of distributions. When selecting the most favourable 
assumption for each of the ten parameters, the probability of 
exceeding the baseline emissions fell to less than 5 %. But 


when combining all least favourable assumptions, this proba¬ 
bility rose to 82 %. Note that both cases are fully realistic as 
the modeller had no a priori knowledge on the variability of 
these parameters. The choice of deciding between the opti¬ 
mistic and pessimistic assumptions is left to the decision¬ 
maker at the interpretation stage. 

We see that the Monte Carlo method and the Confidence 
Index yield very similar results at high levels of probability. 
This is primarily related to the fact that the same distribu¬ 
tions were selected for the pdfs and the fuzzy intervals. 
However, what we see with the IRS result is the conse¬ 
quence of the incomplete character of information 
pertaining to certain parameters. This is seen in the distance 
between the upper and lower probability bounds. In a real 
situation, an excessive distance between these two bounds 
might motivate the decision-maker to increase the informa¬ 
tion base regarding certain critical parameters, in order to 
reduce the uncertainty. Such a decision could not ensue from 
a purely probabilistic calculation based on subjective distri¬ 
butions (despite lack of information), because there is no 
way of distinguishing, in the variability of the calculated 
result, what comes from true randomness and what comes 
from incomplete information. Considering the considerable 
sources of uncertainty in LCA, it is felt that it would be 
more faithful to convey, in addition to an indicator for 
decision-making, an appreciation of the extent of the knowl¬ 
edge gaps and their consequences. 

This study used a rather simple case study with only one 
impact category as the focus was put on the methodology. It 
should be noted that the exact same methodology can be 
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applied at the characterisation, normalisation and weighting 
steps, e.g. to propagate uncertainties in characterisation fac¬ 
tors. It is however acknowledged that the implementation of 
such a method in complex systems and in LCA software 
would require substantial computation power. Indeed the 
calculations involve an optimisation step over several pa¬ 
rameters at each run. 


6 Conclusions 

This paper underlines the difference between different types 
of uncertainty in the context of LCA modelling and illus¬ 
trates a methodology that allows such uncertainties to be 
propagated through the LCA model. Rather than to arbi¬ 
trarily select a given mode of uncertainty representation, it is 
proposed that the investigator first considers the information 
that is available and then selects the formalism that seems 
best suited to convey this information. This sets the focus on 
available infonnation and the importance of gathering infor¬ 
mation that is both reliable and technically feasible, rather 
than disguising imprecise information as precise variability. 
If available infonnation is rich, then a purely statistical 
representation mode is in order, but if it is scarce, then it 
may be better conveyed by possibility distributions. The two 
bounding distributions obtained as a result reflect the in¬ 
complete character of the infonnation pertaining to certain 
parameters: one is the “optimistic” distribution obtained 
when using all favourable values of input possibility distri¬ 
butions, the other one is the “pessimistic” distribution. 
Finally, at the interpretation step, a single distribution can 
be computed by assigning weighs to these two bound dis¬ 
tributions, reflecting the decision-maker’s aversion to risk. 
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Appendix: Calculation of the global warming (GW) 
impact 


Cultivation and harvest of willow (life cycle of 21 years) 

C0 2 _in — (Cem_y2 13 X Cj n _ cu ltivation ^ X C in_harvest) /2 1 

x 44/12 

( 4 ) 

N 2 O em = (N 2 O d + N 2 Oj) x 44/28 x N 2 OCF/l, 000 (5) 


Where: 

C0 2 in Yearly (average) C0 2 emissions from culti¬ 

vation of willow Mg C0 2 ha -1 year -1 


cultivation 
Cin harvest 


V'cm yr 2 

N 2 O em 

N 2 O d 

N 2 Oi 

N 2 OCF 


Yearly net uptake of carbon during the 13 
cultivation years Mg C ha 1 year 1 
Yearly net uptake of carbon during the 5 
harvest years (this parameter being strongly 
correlated to C in cultivation it is later replaced 
by c in -cultivation - 0.78) Mg C ha -1 year -1 
Emissions of carbon during year 2 (assumed 
equal to 5.32) Mg C ha -1 year -1 
Yearly emissions of N 2 0 from cultivation of 
willow Mg C0 2 -eq ha -1 year -1 
Yearly direct emissions of N 2 0 from 
cultivation of willow Mg N ha 1 year 1 
Yearly indirect emissions of N 2 0 from 
cultivation of willow Mg N ha -1 year -1 
Characterisation factor of N 2 0 for GW 
kg C0 2 -eq/kg N 2 0 


Cultivation and harvest of barley 
C0 2 _ b = -C in _b x 44/12 + Y b x C b x 44/12 (6) 


N 2 O b = (N 2 O d b + N 2 Oi_ b ) x 44/28 

xN 2 OCF/l,000 (7) 

Where: 

C0 2 b Yearly C0 2 emissions from cultivation and 
harvest of barley Mg C0 2 ha 1 year 1 
Cm b Yearly net uptake of carbon during cultivation 
and harvest of barley Mg C ha 1 year 1 
7b Yield of cultivation of barley (at the field gate) 

Mg DM ha -1 year -1 

Cb Carbon content of barley %DM 

N 2 Ob Yearly emissions of N 2 0 from cultivation and 
harvest of barley Mg C0 2 -eq ha -1 year -1 
N 2 O d b Yearly direct emissions of N 2 0 from cultivation 
and harvest of barley Mg N ha 1 year 1 
N 2 0; b Yearly indirect emissions of N 2 0 from cultivation 
and harvest of barley Mg N ha 1 year -1 
N 2 OCF Characterisation factor of N 2 0 for GW 
kg C0 2 -eq/kg N 2 0 

Co-firing 

CF = Yield x C w x 44/12 (8) 

Where: 

CF Yearly C0 2 emissions from co-firing of willow 
Mg C0 2 ha -1 year -1 

Yield Yield of cultivation of willow (at the field gate) 

Mg DM ha -1 year -1 

C w Carbon content of willow %DM 
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Avoided energy production 

EP = Yield x (LHV x (1 — Loss) — watercontent/(l — watercontent) 
xtEOheating) x (elec_rec/3.6 x GHG e i e c + heat_rec x GHGheat) 


(9) 


Where: 


EP 

Yearly avoided GHG emission from energy 
production Mg C0 2 -eq ha~' year -1 

Yield 

Yield of cultivation of willow (at the field 
gate) Mg DM ha~' year -1 

LHV 

Lower heating value of willow as dry matter 
GJ Mg 1 DM 

Loss 

Loss of carbon during drying and storage of 
willow % 

Watercontent 

Water content of willow after field drying % 

H 2 Oheating 

Energy needed for water content 
evaporation GJ Mg -1 

elec_rec 

Electricity recovery from LHV % 

heat rec 

Heat recovery from LHV % 

GHG elec 

GHG emissions from electricity production 
in DK Mg C0 2 -eq MWL 1 

GHG heat 

GHG emissions from heat production in DK 
Mg C0 2 -eq GJ -1 


Total net impact 

TNI = iLUC + 20 

x(C0 2 _ m + N 2 O em - (C0 2 _ b + N 2 O b ) + CF - EP) 

( 7 ) 

Where: 

TNI Total net impact on GW over 20 years 
Mg C0 2 -eq ha -1 year -1 

ILUC Indirect land use change Mg C0 2 -eq ha -1 year 1 

C0 2 i n Yearly C0 2 savings from cultivation and harvest 
of willow (average) Mg C0 2 ha 1 year 1 

N 2 O em Yearly emissions of N 2 0 for willow cultivation 
Mg C0 2 -eq ha -1 year -1 

C0 2 b Yearly C0 2 savings from cultivation and harvest 
of barley Mg C0 2 ha -1 year -1 

N 2 O h Yearly emissions of N 2 0 for barley cultivation 
Mg C0 2 -eq ha -1 year -1 

CF Yearly C0 2 emissions from co-firing of willow 
Mg C0 2 ha 1 year 1 

EP Yearly avoided GHG emission from energy 
production Mg C0 2 -eq ha~' year -1 


References 

Andre JCS, Lopes DR (2012) On the use of possibility theory in 
uncertainty analysis of life cycle inventory. Int J Life Cycle 
Assess 17:350-361 

Ardente F, Beccali M, Cellura M (2004) F.A.L.C.A.D.E.: a fuzzy 
software for the energy and environmental balances of products. 
Ecol Model 176:359-379 

Baudrit C, Guyonnet D, Dubois D (2005) Post-processing the hybrid 
approach for addressing uncertainty in risk assessments. Environ 
Eng 131:1750-1754 

Baudrit C, Dubois D, Guyonnet D (2006) Joint propagation and 
exploitation of probabilistic and possibilistic information in risk 
assessment models. IEEE Trans Fuzzy Syst 14:593-608 
Benetto E, Dujet C, Rousseaux P (2008) Integrating fuzzy multicriteria 
analysis and uncertainty evaluation in life cycle assessment. 
Environ Model Softw 23:1461-1467 
Chevalier J-L, Le Teno JF (1996) Life cycle analysis with ill-defined 
data and its application to building products. Int J Life Cycle 
Assess 1:90-96 

Clavreul J, Guyonnet D, Christensen TH (2012) Quantifying uncer¬ 
tainty in LCA-modelling of waste management systems. Waste 
Manage 32:2482-2495 

Couso I, Moral S, Walley P (2000) A survey of concepts of indepen¬ 
dence for imprecise probabilities. Risk Decis Policy 5:165-181 
Cruze N, Goel PK, Bakshi BR (2013) On the “rigorous proof of fuzzy 
error propagation with matrix-based LCI”. Int J Life Cycle Assess 
18:516-519 

Dalgaard R, Schmidt JH, Halberg N, Christensen P, Thrane M, Pengue 
WA (2008) LCA of soybean meal. Int J Life Cycle Assess 
13:240-254 

Danish Energy Agency, energinet.dk (2010) Technology Data for 
Energy Plants. Danish Energy Agency, Copenhagen, Denmark. 
http://www.ens.dk/Documents/Netboghandel%20- 
%20publikationer/2010/Technology _data_for_energy_plants.pdf 
. Accessed 4 December 2012 

DONG Energy A/S, Energinet.dk, Vattenfall A/S (2010) 
Livscyklusvurdering - Dansk el og kraftvarme (In Danish). 
http://energinet.dk/DA/KLIMA-OG-MlLJOE/ 
Livscyklusvurdering/Sider/LCA-resultater-for-dansk-el-og- 
kraftvarme-2008.aspx. Accessed 3 December 2012 
Dubois D (2006) Possibility theory and statistical reasoning. Comput 
Stat Data Anal 51:47-69 

Dubois D, Prade H (1988) Possibility theory. Plenum, New York 
Dubois D, Prade H (2008) Possibility theory: an approach to comput¬ 
erized processing of uncertainty. Plenum, New York 
Dubois D, Prade H (2009) Formal representations of uncertainty. In: 
Bouyssou D, Dubois D, Pirlot M, Prade H (eds) Decision-making 
process-concepts and methods. Chapter 3. London: ISTE & 
Wiley, pp 85-156 

Dubois D, Guyonnet D (2011) Risk-infonned decision-making in the 
presence of epistemic uncertainty. Int J Gen Syst 40:145-167 
Edwards R, Mulligan D, Marelli L (2010) Indirect land use change 
from increased biofuels demand. Comparison of models and re- 


Springer 




Int J Life Cycle Assess (2013) 18:1393-1403 


1403 


suits for marginal biofuels production from different feedstocks. 
Luxembourg: Publications Office of the European Union, http:// 
ec. europa.eu/energy /renewables/studies/doc/1 and_use_change/ 
study_4_iluc_modelling_comparison.pdf. Accessed 27 
September 2012 

Energistyrelsen (2011) Forudsaetninger for samfundsokonomiske 
analyser pa energiomradet (In Danish). Danish Energy Agency, 
Copenhagen, Denmark. http://www.ens.dk/da-DK/Info/ 
TalOgKort/Fremskrivninger/beregningsforudsatninger/ 
Documents/Foruds%C3%A6tninger%20for%20 
samfunds%C3%B8konomiske%20analyser%20p%C3%A5%20 
energiomr%C3%A5det%202011.pdf. Accessed 27 September 
2012 

Ferson S, Ginzburg LR (1996) Different methods are needed to prop¬ 
agate ignorance and variability. Reliability Eng Syst Saf 54:133- 
144 

Frischknecht R, Jungbluth N, Althaus HJ, Doka G, Dones R, Heck T, 
Hellweg S, Hischier R, Nemecek T, Rebitzer G, Spielmann M 
(2005) The ecoinvent database: overview and methodological 
framework. Int J Life Cycle Assess 10:3-9 
Gonzalez B, Adenso-Diaz B, Gonzalez-Torre PL (2002) A fuzzy logic 
approach for the impact assessment in LCA. Resour Conserv 
Recycl 37:61-79 

Guereca LP, Agell N, Gasso S, Baldasano JM (2007) Fuzzy approach 
to life cycle impact assessment—an application for biowaste 
management systems. Int J Life Cycle Assess 12:488—496 
Guyonnet D, Bourgine B, Dubois D, Fargier H, Come B, Chiles JP 
(2003) Hybrid approach for addressing uncertainty in risk assess¬ 
ments. Environ Eng 129:68-78 

Hamelin L, Jorgensen U, Petersen BM, Olesen JE, Wenzel H (2012) 
Modelling the carbon and nitrogen balances of direct land use 
changes from energy crops in Denmark: a consequential life cycle 
inventory. Glob Change Biol Bioenergy 4:889-907 
Heijungs R, Tan RR (2010) Rigorous proof of fuzzy error propagation 
with matrix-based LCI. Int J Life Cycle Assess 15:1014-1019 
Hong J, Shaked S, Rosenbaum RK, Jolliet O (2010) Analytical uncer¬ 
tainty propagation in life cycle inventory and impact assessment: 
application to an automobile front panel. Int J Life Cycle Assess 
15:499-510 

Huijbregts MAJ, Gilijamse W, Ragas AMJ, Reijnders L (2003) 
Evaluating uncertainty in environmental life-cycle assessment. A 
case study comparing two insulation options for a Dutch one- 
family dwelling. Environ Sci Technol 37:2600-2608 
Hurwicz L (1951) Optimality criteria for decision making under igno¬ 
rance. Cowles Commission discussion paper, Statistics No. 370 
Imbeault-Tetreault H, Jolliet O, Deschenes L, Rosenbaum RK (2013) 
Analytical propagation of uncertainty in life cycle assessment 
using matrix formulation. J Ind Ecol. doi: 10.1111/jiec. 12001 
Lindley DV (1971) Making decisions. Wiley-Interscience, London 
Lloyd SM, Ries R (2007) Characterising, propagating and analyzing 
uncertainty in life-cycle assessment, a survey of quantitative 
approaches. J Ind Ecol 11:161-179 


Morgan MG, Henrion M (1990) Uncertainty: a guide to dealing with 
uncertainty in quantitative risk and policy analysis. Cambridge 
University Press, New York 

Reap J, Roman F, Duncan S, Bras B (2008) A survey of 
unresolved problems in life cycle assessment, Part 1: goal 
and scope and inventory analysis. Int J Life Cycle Assess 
13:290-300 

Schmidt JH (2008) System delimitation in agricultural consequen¬ 
tial LCA—outline of methodology and illustrative case study 
of wheat in Denmark. Int J Life Cycle Assess 13:350-364 
Searchinger TD (2010) Biofuels and the need for additional carbon. 
Environ Res Lett 5:024007 

Searchinger TD, Heimlich R, Houghton RA, Dong F, Elobeid A, 
Fabiosa J, Tokgoz S, Hayes D, Yu TH (2008) Use of U.S. 
croplands for biofuels increases greenhouse gases through emis¬ 
sions from land-use change. Science 319:1238-1240 
Shafer G (1976) A mathematical theory of evidence. Princeton 
University Press 

Shulman N, Feder M (2004) The uniform distribution as a universal 
prior. IEEE Trans Inform Theory 50:1356-1362 
Sonnemann GW, Schuhmacher M, Castells F (2003) Uncertainty as¬ 
sessment by a Monte Carlo simulation in a life cycle inventory of 
electricity produced by a waste incinerator. J Cleaner Prod 
11:279-292 

Tan R (2008) Using fuzzy numbers to propagate uncertainty in matrix- 
based LCI. Int J Life Cycle Assess 13:585-592 
Tan R, Culaba AB, Purvis MRI (2004) POLCAGE 1.0—a possibilistic 
life-cycle assessment model for evaluating alternative transporta¬ 
tion fuels. Environ Model Softw 19:907-918 
Thabrew L, Lloyd S, Cypcar CC, Hamilton JD, Ries R (2008) Life 
cycle assessment of water-based acrylic floor finish maintenance 
programs. Int J Life Cycle Assess 13:65-74 
Tonini D, Astrup T (2012) Life-cycle assessment of biomass-based 
energy systems: a case study for Denmark. Appl Energy 99:234— 
246 

Tonini D, Hamelin L, Wenzel H, Astrup T (2012) Bioenergy produc¬ 
tion from perennial energy crops: a consequential LCA of 12 
bioenergy scenarios including land use changes. Environ Sci 
Technol 46(24)43521-13530 

Weckenmann A, Schwan A (2001) Environmental life cycle assess¬ 
ment with support of fuzzy-sets. Int J Life Cycle Assess 6:13-18 
Weidema B (2003) Market information in life cycle assessment. 
Environmental Project No. 863 2003 Miljoprojekt. Danish 
Environmental Protection Agency, http://www2.mst.dk/udgiv/ 
publications/2 003/87-7972-99 l-6/pdf/8 7-7 972-992-4.pdf. 
Accessed 27 September 2012 

Weidema B, Frees N, Nielsen AM (1999) Marginal production tech¬ 
nologies for life cycle inventories. Int J Life Cycle Assess 4:48- 
56 

Williams E, Weber C, Hawkins T (2009) Hybrid approach to 
managing uncertainty in life cycle inventories. J Ind Ecol 
15:928-944 


^) Springer 



